Kondo destruction and valence fluctuations in an Anderson model 



(N 

O 

(N 

o 



J. H. Pixlcy,^ Stefan Kirchner,^'^ Kevin Ingersent,"' and Qimiao Si^ 

^Department of Physics & Astronomy, Rice University, Houston, Texas, 77005, USA 
^Max Planck Institute for the Physics of Complex Systems, 01187 Dresden, Germany 

''Max Planck Institute for Chemical Physics of Solids, 01187 Dresden, Germany 
^Department of Physics, University of Florida, Gainesville, Florida 32611-8440, USA 

(Dated: July 11, 2012) 

Unconventional quantum criticality in heavy-fermion systems has been extensively analyzed in 
terms of a critical destruction of the Kondo effect. Motivated by a recent demonstration of quantum 
criticality in a mixed- valent heavy-fermion system, /3-YbAlB4, we study a particle-hole-asymmetric 
Anderson impurity model with a pseudogapped density of states. We demonstrate Kondo destruc- 
tion at a mixed- valent quantum critical point, where a collapsing Kondo energy scale is accompanied 
by a singular charge-fluctuation spectrum. Both spin and charge responses scale with energy over 
temperature (i^/T) and magnetic field over temperature (H/T). Implications for unconventional 
quantum criticality in mixed- valence heavy fermions are discussed. 

PACS numbers: 71.10.Hf, 71.27.-|-a, 75.20.Hr 
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Competing interactions in quantum systems give rise 
to zero-temperature phase transitions. If it is continuous, 
such a transition takes place at a quantum critical point 
(QCP). There is mounting evidence, especially in heavy- 
fermion systems, that a QCP can underlie unconven- 
tional superconductivity related considerations have 
been applied to high-temperature cuprate and iron pnic- 
tidc superconductors 0- It is standard to describe a 
QCP within the Ginzburg-Landau- Wilson (GLW) frame- 
work: critical destruction of an order parameter charac- 
terizing a spontaneously broken symmetry gives rise to 
collective modes associated with order-parameter fluctu- 
ations In the context of antiferromagnetic metals, 
this is referred to as a spin-density- wave QCP 

Recent experiments in heavy-fermion metals have 
clearly established the existence of a novel class of anti- 
ferromagnetic QCPs, characterized by non-Fermi liquid 
behavior and a; /T scaling in the dynamical spin suscepti- 
bility . There are indications that such unconventional 
QCPs also promote superconductivity 0. These QCPs 
defy description in terms of a GLW functional [tI, 
their understanding requires the introduction of quan- 
tum modes beyond order-parameter fluctuations. The 
proposed additional modes arc associated with the criti- 
cal destruction of the Kondo effect In the paramag- 
netic phase, Kondo singlets form and generate Kondo res- 
onances, thereby turning the local moments into single- 
electronic excitations and enlarging the Fermi surface. 
The destruction of the Kondo effect across the antiferro- 
magnetic QCP suppresses the Kondo resonances, mak- 
ing the Fermi surface small. Critical Kondo destruction 
therefore manifests itself in a discontinuous evolution of 
the Fermi surface across the transition, as has been ob- 
served through quantum oscillation and Hall effect mea- 
surements [J^]. 

Theoretical studies of critical Kondo destruction have 
largely been confined to the Kondo-lattice limit of inte- 



ger valence. In rare-earth intermetallics, superconduc- 
tivity is believed also to arise in the vicinity of valence 
transitions which have been found to be first or- 
der. Until recently, there has been no significant evi- 
dence for a QCP associated with valence fluctuations. 
The situation has changed with the observation of mixed 
valency in the ytterbium-based heavy-fermion supercon- 
ductor /3-YbAlB4 lull, which is quantum critical under 
ambient conditions jl2l |. In an applied magnetic field, the 
magnetization obeys H/T scaling [13[, consistent with 



the uj/T scaling seen previously near the unconventional 
QCPs of antiferromagnetic heavy-fermion compounds. 
These properties implicate /3-YbAlB4 as a strong can- 
didate for a mixed- valent heavy-fermion QCP, and raise 
the prospect that the material's unusual scaling behavior 
can be understood in terms of critical Kondo physics. 

At first glance, critical Kondo destruction at mixed 
valence appears unlikely. Kondo destruction in a Kondo 
lattice amounts to the localization of / electrons. While 
unconventional, this is physically transparent, because 
localization can readily arise for a commensurate filling of 
an electronic orbital (one /-electron per site). At mixed 
valence, the situation is more subtle because the / or- 
bital has a fractional, generally incommensurate, per-site 
occupancy, and there is no mechanism known for elec- 
tron localization at incommensurate fillings. This leads 
to important questions of principle: Can critical Kondo 
destruction occur in the presence of valence fluctuations 
and, if so, how does the criticality compare to its local- 
moment counterpart? For instance, are charge excita- 
tions part of the critical fluctuation spectrum? 

In this Letter, we address these issues in the mixed- 
valence regime of an Anderson impurity model whose 
conduction-electron density of states features a pseudo- 
gap centered on the Fermi energy. We focus on an im- 
purity model because of the local nature of the Kondo- 
destruction physics; formally, the Kondo lattice model 
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FIG. 1: (a) Binder cumulant B{U,P) vs. U for r = 0.6, To = 
O.ID, and ed = —0.05D, and the labeled temperatures. Error 
bars were obtained from a jackknife analysis of the CT-QMC 
data, (b) Blow up of the same data around the intersection 
of curves, which determines Uc/D = 0.063125 ± 0.0008. 



can be treated through an effective impurity model in the 
extended dynamical mean field approach IJ, [15[ . Given 
that the commensurate-filling (i.e., local- moment) limit 
of the model exhibits critical Kondo destruction and asso- 
ciated dynamical scaling properties [lij , we consider the 
pseudogappcd density of states to provide a prototype 
setting to search for a Kondo-destruction QCP at mixed 
valence. Our model has the advantage of being amenable 
to study using reliable methods: the continuous-time 
quantum Monte Carlo (CT-QMC) method 1171 and the 

numerical renormalization group (NRG) [iilillii- 

Surprisingly, we do find critical Kondo destruction in 
this mixcd-valent model. The critical properties in the 
spin sector reflect the collapse of an energy scale as the 
QCP is approached from the Kondo-screened side but not 
from the Kondo-dcstroyed side, much as in the integer- 
valent (local- moment) limit. By contrast, the charge sec- 
tor shows a collapsing energy scale on both sides of the 
QCP. The critical point displays H/T (and w/T) scal- 
ing. This existence proof for a Kondo destruction QCP 
at mixed valence makes it feasible to interpret the H/T 
and related scaling properties of /3-YbAlB4 in terms of an 
interacting fixed point. We note that the same model is 
also relevant to impurity physics in d-wave superconduc- 
tors and graphcne, where the density of bulk fermionic 
states goes to zero at the chemical potential [20| . 

The Anderson impurity Hamiltonian is 



H = '^[ekcl^Cka + V {dicker + ii-c.)] +ednd + Und^ndi, 

(1) 

where Cka annihilates a conduction-band electron of en- 
ergy Cfc, d^r annihilates an electron of energy Ed in the im- 
purity level, U is the electron-electron repulsion within 
the impurity level, V is the hybridization taken to be mo- 
mentum independent, Uda = d^d^, and nd = Ud^ + Udi- 
The band density of states vanishes in a power-law fash- 
ion at the Fermi energy (ep = 0): 

p(e) - ^ -5(6 - ek) = PoWDl^eiD - \e\). (2) 



The impurity-band interaction is completely specified by 



FIG. 2: Valence and local spin properties vs. u = U/Uc — 1 for 
r — 0.6, Fo ~ 0.1-D, and Sd = — 0.05-D: (a) Occupancy (ud) 
at the labeled temperatures [2^. The QCP {u = 0) occurs at 
mixed valence, i.e., (ud) 7^ 1. (b) Local magnetization Mioc, 
showing quenching of the impurity spin for u < but the 
emergence of a free local moment for u > 0. 



the imaginary part of the hybridization function, r(e) ~ 
^ Ek V^H^ - ^h) = ro\e/D\\ where Pq = npoV^. 

The critical properties of the model with particle-hole 
(p-h) symmetry (e^ = — [7/2) and its Kondo limit {U » 
Po, where local charge fluctuations are negligible) have 
been investigated in a number of analytic and numerical 
studies [IH-il. The breaking of p-h symmetry 



is irrelevant for pseudogap exponents r in the range < 
r < r* ~ 0.375, but becomes relevant for r > r*, leading 
to a mixed- valent QCP [l9j : r = 1 serves as an upper 
critical "dimension" , above which the critical properties 
have a mean-field character 1^ 2^ . 

Here, we investigate the p-h-asymmetric pseudogap 
Anderson model by varying U for fixed Pq and Ed to pass 
from a Kondo-screened strong-coupling phase {U < Uc) 
to a Kondo-destroyed local-moment phase {U > Uc). 
We apply the CT-QMC technique, which was recently 
shown to be able to reach temperatures T sufficiently 
low to access the quantum critical regime We 
measure the dynamical local spin and charge suscep- 
tibilitics, x.(t,/3) = {TrS,iT)SM) and Xc(r,/3) = 
(T7- :n(T): :n(0):), respectively, where Sz = ■^{n^ — ni), 
- (Ea'^'y)' and /3 = 1/T (taking ks = 1) 
plays the role of the system size. The corresponding 

static susceptibilities follow from Xc,s(/3) = Jq dTXc,s{T). 
Measuring powers of the local magnetization (Af") = 

([■^ /o^ d'TSz{T)Y') allows construction of the Binder cu- 
mulant [23 -B(C/,/3) = {Mf)l{Mlf. We supplement our 
T > (finite-/3) CT-QMC results with static quantities 
calculated arbitrarily close to T = (/? = 00) using the 
NRG method as adapted to treat pseudogap impurity 
problems 16, 3 l^- NRG results presented below were 
obtained with Wilson discretization parameter A = 9, 
with Po corrected [31 to compensate for the band dis- 
cretization, and retaining all many-body states up to 50 
times the effective bandwidth of each iteration. 

We focus our discussion on the representative case of 
a pseudogap exponent r = 0.6 with Po = 0.1-D and 
£d ~ —0.05-D. Figure [T] plots the variation of the Binder 
cumulant with \J at different temperatures. For small 



3 



U , charge fluctuations are strong and the Binder cu- 
mulant hes above the range 1 < B{U,13) < 3 obeyed 



by a pure spin system [23|, |2J]; in this limit, the low- 



energy behavior is close to that of a pseudogap resonant 
level. As U increases, charge fluctuations are suppressed 
leading at low temperatures to B{U,f3) < 3; this part 
of the strong-coupling phase exhibits a true pseudogap 
Kondo effect. For very large U, by contrast, B{U, (3) 
tends towards 1 at low temperatures, suggesting the pres- 
ence of a decoupled impurity spin, characteristic of the 
local-moment phase. We locate the phase boundary by 
the intersection of B{U, (3) curves for different tempera- 
tures m at Uc/D = 0.06313 ± 0.0008. The NRG gives 
Uc = 0.064501^, a small shift that can likely be attributed 
to residual effects of NRG discretization. The mixcd- 
valent nature of the QCP is demonstrated in Fig. [D^a), 
where the local occupation (ud) is seen to differ from 
unity at U ^ Uc- Note also that (rid) displays signiflcant 
temperature dependence in the vicinity of the QCP. 

We are now in a position to look for a critical destruc- 
tion of the Kondo effect in this mixed- valent QCP, i.e., 
the continuous vanishing of an effective Kondo energy 
scale signaled by the divergence of the zero-tempcrature 
static local spin susceptibility Xs as U approaches Uc 
from below. Such a divergence is indeed seen in our 
zero-temperature Xs vs. U data [Fig. [3l[a)] and in the 
temperature dependence of Xs at ?7 = Uc [Fig. Eljb)]. 
Figure [lib) shows the U dependence of the local mag- 
netization Mioc = limft_j.o limT-i.o(-^'^z), where /i is a lo- 
cal magnetic field entering a term h{n-f — n^)/2 (with 
gfiB = 1) added to Eq. ([T|). Since Mioc = throughout 
the strong-coupling phase, and Mioc rises continuously 
from zero on entry to the local-moment phase, this quan- 
tity serves as an order parameter for the quantum phase 
transition. Our results can be summarized as 



XsiT,U^Uc) 
Xs{T^Q,U<Uc) 
Mioc{T = 0,U>Uc) 



T- 



(3) 



where u = U/Uc - 1. We find Xs = 0.80(3) from CT- 
QMC, in excellent agreement with the NRG value Xs = 
0.7908(3); the NRG also yields 7^ = 1.42(2) and = 
0.1874(2). These power-law behaviors are all defining 
characteristics of critical Kondo destruction. 

To probe valence fluctuations near the QCP, we turn 
to the static local charge susceptibility Xc{T,U). As 
shown in Fig. |4l[a), Xc{T = 0,U) increases with U in the 
strong-coupling phase and diverges as U ^ in a 
manner similar to Xs{T = 0,U). In the local-moment 
phase, the spin and charge responses are very different: 
Xs{T = 0,U) = 00 but Xc{T = 0,U) remains flnite, al- 
though it diverges as U ^ Uc ■ In other words, the va- 
lence fluctuation energy scale is nonzero in both phases, 
vanishing only when U approaches Uc from either side. 
At U = Uc, Xc has a singular temperature dependence as 
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FIG. 3; Local static spin susceptibility Xs{T,U) for r = 0.6, 
To = O.ID, and ea = -0.05L>, (a) vs. u = U/Uc-l at the 
labeled temperatures [2^, and (b) vs. T at the critical point 
U^Uc. 



shown in Fig. IHb). These behaviors are consistent with 



Xc(T, U = Uc) 
Xc{T = 0,U) 



T- 



(4) 



CT-QMC yields Xc = 0.36(3), while the NRG gives Xc = 
0.120(1) (extracted at temperatures much lower than can 
be accessed by CT-QMC) and 7c = 0.21(1). The differ- 
ence between the two Xc values stems from a very slow 
crossover to the quantum critical regime [Fig.m^b) inset]. 
The much wider crossover window for Xc compared with 
Xs [Fig-EIb)] likely arises because Xc < Xg, meaning that 
lower temperatures must be reached before sub-leading 
contributions to Xc become negligible. We stress that the 
singularity of the charge susceptibility is unique to the 
mixed-valence QCP, and does not appear at its integer- 
valence counterpart. 

We reach the important conclusions that the Kondo 
destruction occurs at a genuinely mixed- valent QCP and 
that valence fluctuations are part of the critical spectrum. 
Calculations for level energies £d ^ — 0.05D (and hence 
different critical occupancies {ud)) indicate that the crit- 
ical exponents deflned above depend on the band expo- 
nent r, but not on the impurity valence. This implies 
that the divergence of the static charge susceptibility is 
a universal property. At the same time, we flnd that 
the critical behavior in the spin sector coincides with the 
model in its integer valence limit, i.e., the p-h-asymmetric 
pseudogap Kondo model 

We now discuss the dynamical scaling of Xs{t,T) and 
Xcir, T). In analogy with the spin response at the Kondo 
destruction QCP in the usual Kondo limit [2^, [2^, we 
flnd that at U = Uc both Xs{t,T) and Xc(t, T) collapse 
onto the conformal scaling form, showing a power-law de- 
pendence on ttT/ sin(7rrT) with exponents rjc^s- For the 
temperatures considered, the charge susceptibility has 
not yet reached its asymptotic power-law behavior [based 
on Fig. m^b)]. Our results thus imply that both lead- 
ing and sub- leading terms of the critical Xc{t, T) scale in 
terms of ttT/ sin(7rrT). The sca ling form means Xs(w, T) 
and Xc{'^, T) obey uj/T scaling [23| at U = Uc- 

We next consider the effect on the QCP of applying 
a finite local magnetic field h. Consistent with the lo/T 
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FIG. 4: Local static charge susceptibility Xc{T, U) for r = 
0.6, To = Q.ID, and = -0.05-D, (a) vs. u = (7/C/c-l at 
the labeled temperatures [2^, and (b) vs. T at the critical 
point U — Uc, where the discrepancy between CT-QMC and 
NRG data is due mainly to the difTerence in Uc values. Inset; 
Xc{T, Uc) over a wider range of T, showing the slow crossover 
behavior. 
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5: Scaling of critical spin and charge responses for 
r = 0.6, Fo = 0.1-D, Sd = -0.05D and U ~ Uc. (a,b) 
Dynamical spin and charge susceptibilities vs. nT/ sm{nTT). 
Both susceptibilities show excellent scaling collapse over two 
decades of 7rT/sin(7rTT), with ris = 0.20(3) and r?^ = 0.67(3). 
(c,d) Static spin and charge susceptibilities vs. T/h in the 
labeled local magnetic fields h. 



scaling we find T/h scaling for fields \h\ < Tk, i.e., 

XcAT,h,Uc)^h-y^-Us{T/h), 



(5) 



where fc,s{x) is a scaling function [see Figs. [SJc) and 
[5jd)]. Scaling collapse in T/h further reflects the in- 
teracting nature of this mixed-valence QCP. At such an 
interacting QCP, critical exponents satisfy hyperscaling 
relations that imply j/c,s = a;c,s, equalities confirmed by 
our results to within numerical accuracy. 

Our results are to be contrasted with the generaliza- 
tion__of the spin-density-wave QCP to_the valence sec- 
such a mixed- 



tor 



271 128[ . Like its spin counterpart 
valent QCP is Gaussian (noninteracting) and is not ex- 
pected to obey either energy-over-temperature or field- 
over-temperature scaling. 

We now briefly consider the case r = 1, motivated 
by critical Kondo screening in graphene and d-wave su- 



perconductors. For Fq = 0.1-D and Sd = — 0.05-D, we 
find Uc/D = 0.05475 ±0.0006 with CT-QMC, or 0.05562 
with NRG. Both methods indicate that the local static 
charge susceptibility diverges at along with the spin 
susceptibility. We find Xc to be an increasing function 
of r, but logarithmic corrections to scaling (l6[ prevent 
reliable determination of critical exponents for r = 1. 

This work provides new insights into the unusual crit- 
ical properties of /3-YbAlB4 , suggesting that Kondo 
destruction can occur in this material even though it is 
mixed valent. (Mixed valence is natural in this material 
given that its onset Kondo temperature is high — on the 
order of 200 K — and its mass enhancement is moderate.) 
The demonstration of h/T scaling provides evidence that 
the experimentally observed field-over-temperature scal- 
ing signals Kondo destruction. Our finding of a rapid 
variation of ( nd) near the QCP, which can be tested ex- 
perimentally [29[, suggests that the concentrated lattice 
system is essentially quantum critical over a range of den- 
sities, leading to the exciting possibility of quantum crit- 
icality occurring over a region of parameter space rather 
than just at an isolated point. Finally, our work raises 
intriguing questions about the extent to which quantum- 
critical magnetic and valence degrees of freedom influence 
the superconductivity observed in /3-YbAlB4. 

From a general theoretical perspective, how quantum 
criticality can go beyond the GLW framework of order- 
parameter fluctuations is a fundamental problem that is 
important not only for heavy-fermion metals but also for 
QCPs arising in insulating magnets and other strongly 
correlated systems. At present, there are few concrete 
theoretical examples for such unconventional QCPs. By 
identifying a new QCP in this category, our results pro- 
vide another setting to gain intuition about bcyond-GLW 
QCPs in general. 

In summary, we have shown that mixcd-valcnt quan- 
tum criticality can display the phenomenon of Kondo 
destruction. The quantum critical point has a collaps- 
ing Kondo energy scale and a singular charge-fluctuation 
spectrum. The valence varies strongly with tempera- 
ture near the critical point. In the concentrated lattice 
case, similar quantum critical behavior is expected to oc- 
cur over an extended range of parameters. Our results 
raise the prospect of unconventional quantum criticality 
in mixed- valent systems beyond /3-YbAlB4. 
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